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ABSTRACT 




Two important capabilities have been added to the spinning-body 
version of the EDDYBL computer program currently operational at the 
NASA Ames Research Center. First, the conventional mixing-length 
model, specialized for thick boundary layers, has been added to 
the program's array of possible turbulence models. Second, pro- 
vision has been added for using a more general model for pressure- 
strain correlation terms in the Reynolds-stress-model option. 
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NOTATION 


Symbol 

A 

<=f 

c 

°1» C 2> 
°3 

°P 

D U 

e 

h 

H(3) 

J 

P U 

P 


Definition 

Van Driest damping coefficient. Equation (42) 

Skin friction 

Constant in the law of the wall. Equation (30) 

Closure coefficients 

Specific heat at constant pressure 

Production rate tensor, Equation (6) 

Turbulent kinetic energy, - x ii 
Specific mean enthalpy 
Heaviside stepfunction 

Absolute value of mean strain-rate tensor 
Production rate tensor. Equation (6) 

Half the trace of production rate tensors. Equation (7) 




9i 

q 

r 


R e ,R 

Rea 

0 k 

s 

Sf> S 


0) 


t 



Laminar, turbulent Prandtl numbers 
Turbulent heat flux vector 
Component of q i normal to surface 
Radial coordinate 

Closure coefficients in viscous modifications 
Reynolds number based on 0 k 
Arc length 

Value of s at the end of, beginning of' transition 

Mean strain-rate tensor 

Time 


iv 


Symbol 

T 

T t 

u i 

u,v,w 

u + 

u 


u e 

V 

X 1 

x,y,z 


Definition 

Mean temperature 
Mean total temperature 
Mean velocity vector 
x,y,z components of u A 

Dimensionless sublayer-scaled velocity, u/u i 
Friction velocity, /r^7pT 


Velocity at boundary-layer edge 

Dimensionless transformed vertical velocity 
Position vector 

Streamwise, normal, lateral components of x ^ 
Dimensionless sublayer-scaled normal distance, u T y/v 
Eddy-viscosity coefficient. Equation (44) 
Launder-Reece-Rodi closure coefficients 
Closure coefficients in dissipation terms 


a 

AAA 

«, e, y 


•V 

B 

y,y* 

y 

r(s) 

6 



e 

e + 

e i ,e o 


Eddy-viscosity coefficient, Equation (49) 

Closure coefficients controlling production 
lntermittancy in boundary layer. Equation (45) 
Intermittancy in streamwise direction, Equation (51) 
Boundary layer thickness 

Kinetic displacement thickness. Equation (46) 
Kronecker delta 

Width of transition region 

Eddy diffusivity 

Dimensionless eddy diffusivity, e/v 
Inner, outer layer eddy diffusivity 


v 


Symbol Definition 

0 k Kinetic momentum thickness, Equation (47) 

•c Kantian' s constant 

X Closure coefficient In viscous modifications 

A 

X Intermittancy coefficient } Equation (52) 

M Molecular viscosity 

v Kinematic molecular viscosity, p/p 

ir Eddy viscosity parameter, Equation (48) 

p Mean density 

cr,o*,o## Turbulent "Prandtl" numbers 

T iJ Reynolds stress tensor 

t Reynolds shear stress, <-u'v'> 

X Strain-rate parameter, ^3 ~ S mn /e w 2 

w Turbulent dissipation rate 

A Turbulent length scale; also mixing length 

A A Mixing length In inner region 

Subscripts 

e Boundary-layer edge 

o, w Surface 

Other Notation 

For a given variable ip: 

*ip>, Long-time-mass-averaged value of 

ij> f Fluctuating part of ij> 
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RECENT IMPROVEMENTS TO 
THE SPINNING BODY VERSION 
OF THE "EDDYBL" COMPUTER PROGRAM 

1. INTRODUCTION 

In order to analyze effects of three-dimensionality on turbulent 
boundary layers, reliable' experimental data are needed. In order 
to help provide such data, the NASA Ames Research Center has been 
developing an experimental facility for measuring properties of 
swirling axisymmetrlc boundary layers . While such boundary lay- 
ers are not fully three dimensional, they at once (a) display 
coupling amongst all six Reynolds stresses similar to the 3-D 
situation and (b) can be computed with 2-D solution methods. 

Concurrent with development of the experimental facility, the NASA 

l-T 

Ames Research Center has sponsored development of a swirling 

h 

version of DCW Industries' EDDYBL computer program. Prior to this 
research project, the program contained the two equation and Reynolds 
stress models devised by Wilcox and Rube s in. ^ 

Since the time of the program's original development, developments 
in turbulence modeling have spurred Interest in expanding the 
capabilities of the spinning version of EDDYBL, viz: 

1. The pressure-strain-correlation closure 

c 

approximations of Launder, Reece and Rod! 
have been found to be more physically sound 
than those used by Wilcox and Rubesln; 
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2. General overall success with refined 

versions of the mixing-length model such 
as that proposed by Aguilar^ provides 
Justification for Including mixing-length 
predictions In any meaningful boundary- 
layer study. 

The purpose of this project has been to Include the Launder-Reece- 
Rodi (LRR) pressure-strain-correlation terms and the Aguilar 
mixing-length model in the program’s cadre of turbulence models 
and closure approximations. 

Section 2 presents details of analysis based on the LRR pressure- 
strain-correlation terms while Section 3 focuses on the mixing- 
length model. Section ^ summarizes results of this study. The 
Appendix presents modifications to code input and output. 


2. REVISED PRESSURE-STRAIN-CORRELATION 
CLOSURE APPROXIMATIONS 

To enhance the generality of the spinning-body version of the 
EDDYBL computer program (hereafter referred to as EDDYBL for 
brevity), one objective of this project has been to include pro- 
vision for using the launder, Reece and Rodi 5 (LRR) pressure-strain- 
correlation closure approximations. Several steps are involved in 
accomplishing this objective, viz: 

1. Stating the complete model; 

2. Establishing values of the closure coefficients; 

3„ Developing viscous modifications; and 

Incorporating revisions in the program. 

Complete details of Steps 1-3 are given in this section. As for 
incorporating the modifications (Step 4), the most significant 
point is the manner in which input and output are affected. We 
now have several new input variables and certain default values 
have been changed. The Appendix summarizes revised input and 
output . 

2.1 THE MODEL EQUATIONS 

For a compressible fluid of density p moving with mass-averaged 
velocity u^, the Reynolds-stress tensor, , and turbulent heat- 
flux vector, q^, are computed from the following equations: 
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where h is mass averaged enthalpy, |J Is molecular viscosity, 

Pr^ is laminar Prandtl number, t is time, xj. is position vector 
and e is eddy diffusivity defined as the ratio of turbulent 
energy, e » - tii/2, to turbulent dissipation rate, w , i.e. , 

c »• e/u (4) 

Similarly, H is the turbulent length scale defined in terms of 
e and co as 

SL ■ e 1//2 /a) 

4 


(5) 


The tensors ^ , Dy and are defined as 
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3u, 
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( 6 ) 


while the quantity P Is half the trace of P.y (and of D.^), viz 

p * I p kk * ? D kk ■ t mn s nm (7) 


A A A 

Finally, the quantities a, 3, g, 0*, 3**> y, y, a, a* , a* 11 
and are closure coefficients whose values are established in 
the following section. 


2.2 VALUES OF THE CLOSURE COEFFICIENTS 

A A A 

2.2.1 Setting the Values of a. 3. and y. 

As shown by Launder, Reece and Rodi, symmetry properties of the 

A A A 

exact pressure-straln-correlation term imply that a, 3 and y 
can be expressed in terms of a single unknown coefficient, C 2 , 
as follows: . 

“ - rr <° 2 + 8 > J 

e * jx <eo 2 - 2 ) ( (8)' 

y - jg (15C 2 - 1) \ 
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In the LRR model, a value Cg * 0,4 appears most appropriate 
which corresponds to 

5 - 0,76, $ ■ 0.11, y * 0.36} WR» Model 

(9) 

Note that we denote the new model as the WR* model. In the WR 

I 

model there Is no single value for Cg so that some of the sym- 
metry properties of the pressure-straln-correlatlon tensor are 

A 

violated. The values for a, etc., are 

a - 3 - 1/2, y * 4/3 i WR Model (10) 


2,2.2 Setting the Value of . 


Turning to the coefficient C^, study of the decay of anisotropic 
turb ilence and Its asymptotic return to isotropy implies a 


value for which, 
rate parameter x * 


i n the WR model , depends upon the strain 
v S nm S mn /$*w 2 according to 


C 1 « 4.5 - 2.5 exp ( — 5X ) } WR Model (11) 


By contrast, C 1 is constant and equal to 1.5 in the LRR model. 
We assume this value applies for the WR* model. 

C ± - 1.5 } WF?* Model (12) 
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2.2.3 Setting the Values of a. q* and o**. 


For the turbulent Prandtl numbers a, a* t a** we elect to use 
the same values as those In the WR model,, wherefore we have 

o « a* « 1/2, - 2; Both Models (13) 


2.2.4 Setting the Value of g. 

Accurate simulation of the decay of homogeneous, Isotropic 
turbulence means we must require 3 and to stand in the 
following ratio? 

$/3* - 5/3 ; Both Models (14) 

2.2.5 Setting the Values of 8* and y. 


Analysis of the incompressible wall layer shows that the 
coefficient y must satisfy 


e/e 


* _ 


2«k ! 


where y* is the (constant) ratio of shear stress, t 
e3u/3y, i.e.. 


X = v # — 

T w 3 y 


and is given by the following equation: 


(15) 

to 


( 16 ) 




l-a)(C 1 -l+a- 23 )-e(C 1 + 2 - 2 i+ 3 )] 


Also, in the wall layer. 


t/e 


(17) 


( 18 ) 
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while the normal stresses are given by 
<u ,2 >/e * ^ |c^ + 2-2a + ej 

>/e ■ 4 [°i - i+ • -A 

>/e ■ 4 [°i - i+ « + »] 


f 2 


t 2 


( 19 ) 


If, as in the WR model, we select 
0 * * 9/100 


( 20 ) 


then it follows that 

1.097 

• 

J 

WR Model 

Y* s 

1.413 

• 

> 

WR* Model 

and y is given by 

13/11 

• 

1 

WR Model 

Y £ 

4/3 

• 

9 

WR* Model 


( 21 ) 


( 22 ) 


Note that in the wall layer we 


t/e •> 


0.31 

0.36 


now have : 

; WR Model 
; WR* Model 


and 

<u f2 >: <v ,2 > : <w* 2 > = 


4; 2 : 3 J WR Model 
4: 2 : 2.6 j WR* Model 


(23) 


( 24 ) 


2.2.6 Setting the Value of 3»». 


Finally, to specify $#* we note that for an Incompressible 
boundary layer the radial heat flux equation simplifies, in 
the wall layer, to a balance between production and dissipation, 

1 . e . , 

9C T 

0 * - p<v ,2 > -- gP- - $##pwq (25) 


so that q follows a gradient diffusion law defined by 


q 



«v«> !M 

w 3 y 


( 26 ) 


Now, since Equation (16) tells us the eddy diffusivity in the 
momentum equation is e = y#e/u>, the classical analogy between 
momentum and heat transfer encourages us to write 


q 


e 

Pr, 


3 M 


T 


a y 


(27) 


which implies [upon Inspection of Equations (16), (19)> (26) and 
(27)3 that 8#* and Pr,p are related by the following equation: 

2Pr. 


8»* 


3C 


I?* [°i- 


1+0-23 


\ 


( 28 ) 


Hence, for the two models we have (demanding that Pr T = 8/9 in 
the wall layer) : 


B*» = 


0.36 

0.29 


WR Model 

(29) 

WR* Model 



2.3 VISCOUS MODIFICATIONS 


Just as with the WR model, we must devise viscous modifications 
for the new model. Integration through the sublayer with no 
viscous modifications shows that, in terras of standard sublayer 
coordinates, the new model yields 

u + *► 7 In y + + C (-30) 

where C * 8.6 for a perfectly smooth wall. As with the WR model 
we introduce the following viscous modifications: 

0 X » | jl - (l-X 2 ) exp (- £+ /R e )l _1 (31) 

c x ■ [i - d-* 2 ) ex P (- E+/R u)] (32) 

where e + ■ pe/ioy , subscript *> denotes value for high Reynolds 
number turbulence, and R e , R^, X are coefficients to be determined. 


The value of X is set by demanding that turbulence dissipation 
exceed turbulence production below the minimum critical Reynolds 
number for the Blasius Boundary layer. This condition is satis- 
fied provided 


X* s .00*50^ 
3y+Hi-a-S) 


(33) 


Hence, we obtain 


X 


1/14 ; WR Model 

2/31 J WR* Model 


(34) 
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To set the values of R e and R w we proceed as in previous analyses 
noting that there is a locus of pairs of values (R e , R M ) which 

W w 

yield C ■ 5.5 for a perfectly smooth wall (Figure 1). The 
optimum pair appears to be 

R e ■ i, R w - 2.67 j WR* Model 

C35) 

which is close to the pair of values used in the WR model, viz, 

R e * 1, R w »3 » WR Model 

(36) 

Figures 2 and 3 compare predicted results for (R . R ) » (1,2.67) 

G UJ 

and (R e , R^) = (0,0.93). Clearly, the larger value of R e is 
preferable. Further increase in R g tends to drive the turbulence 
production away from the experimental data. 
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u/u T 



Figure 2. Velocity profiles in the sublayer. 
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3 . MIXING-LENGTH MODEL 

The second objective of this project has been to Include the 
7 

Aguilar mixing-length model In the program. Three steps are 
Involved In accomplishing this objective, via: 

1. Stating the complete model; 

2. Specifying the transition point and 
width of the transition region; 

3. Incorporating revisions in the program. 

As with the modifications discussed in Section 2, the Appendix 
summarizes changes in program input and output. Complete details 
of Steps 1 and 2 are given in this section which concludes with 
results of a test case. 

3.1 THE MODEL 

This subsection summarizes the Aguilar mixing 
spinning, thick boundary layer. The Reynolds 
turbulent heat-flux vector are written as 

T u ■ 2e ( s u - 5 u) 

x k 

e MV) 
q i Pr T 

where e is the eddy diffusivity. 

To compute e, the model divides the boundary layer into an inner 
and an outer region. The eddy diffusivity in the inner region. 


■length model for a 
stress tensor and 


(37) 

(38) 
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c ±l is 


e 1 - 


( 39 ) 


where |j| Is the magnitude of the mean strain rate and is the 
inner region mixing length. For flow over an axlsymmetric body 


of radius r Q > these quantities are given by 


(H - f) ! + ( 


mV 

3r j 


*i " « r o ln F 


1 - exp 


-r Q ln(r/r 0 ) 


hr 


( 40 ) 


( 41 ) 


where u denotes streamwise velocity, w is swirl velocity, r is 

radial distance, and k Is Karman’s constant. The quantity A is 

the van Driest damping coefficient defined by 
26 v 


w o 


( 42 ) 


where is the value of J at r * r and v denotes molecular 
0 o 

diffusivity . 


In the outer region, the eddy diffusivity, e , is given by 

e * a U 6, * Y 
o w e k 1 


( 43 ) 


where U is freestream velocity and the quantity a is given by 

v? 


.0168 


, Re 0 > 5000 


° m \‘ 0 l 68 (h^) Re k 

' u+ * > • Re e k < 5ooo w 

Also > In Equation (tn u,. 

by quantity y Is the intermlttaney defined 

* ' Q + 5 -5 (r - r 0 )/«j -i 

_ (If 5) 

In Equations (i| 3 ) and ( « 

and momentum thicknesses defined " * k ^ ^ klnetl ° dis Plaoement 

es defined as follows: 

» ^ / 


W 

k J* ^ ~ U/U e ^ dr 

o 


■/' 


U e ^ - u/ M dr 


V 

Finally, th e quanfcit-v «, 4 

’ • “ *» - ,o™, u , 

"here H(g) ls th „ 

by iSWe StePf " n ^lo„ and ( * is deflned 

8 = Re 0 /4 25 - i 

Ir 


(46) 


(47) 


(48) 


(49) 


The two regions are j 0 in^ k 

viscosity be contlnuoua acro ;;; h e ; equi ~ «»• the eddy 
insured by requiring °undary layer, which i s 

e * min (e e ) 

(50) 
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3.2 TRANSITION POINT AND WIDTH OP TRANSITION REGION 
To complete specification of the mixing-length model we must 
establish the location of transition, the width of the transition 
region, and the streamwlse intermit tancy on the transition region. 
This must be done because, unlike the two-equation and RSE models, 
the mixing-length model has no natural way of predicting 

t 

transition location or width. 

8 

Following Dhawan and Naraslmha, we introduce a streamwlse inter- 
mittancy, f(s), defined as follows: 



where s Is streamwlse distance, s fc is the value of s at the 

A 

transition point, and A is defined by 

* ■ < s f - ' <52) 

with s f denoting the value of s at the end of transition (at 
which point r- 0.98). 


Finally, to determine the width of the transition region we can 
either specify it empirically or we can use the correlation 
devised by Dhawan and Naraslmha, viz, 


where 



5Re 


o,« 


As t - s f - s t 


(53) 

(5*0 
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This completes specification of the mixing-length model. In the 
following subsection we present a simple test case. 

3.3 SAMPLE COMPUTATION 

As a preliminary test case, we present results for an Incompressible 
(M e ■ 0.13) boundary layer on a cylinder of sufficient radius for 
the boundary layer to be essentially two-dimensional (5/r* 1/10) . 
Figure 4 compares skin friction, c f , (a) computed with the mixing- 
length model, (b) computed with the Wilcox-Rubesln RSE model and 

Q 

(c) as correlated by von Karman. As shown, except very near 
transition, the computations and the correlation are within 10# . 

The only significant observation about the program's operation Is 
that, in order to achieve an accurate solution with the mixing- 
length model, the solution should be required to differ from 
iteration to iteration by a tenth of a percent as opposed to one 
percent for the more advanced models. 


MIXING-LENGTH MODEL 
WILCOX-RUBESIN RSE MODEL 
EXACT , BLASIUS FORMULA 

VON KARMAN CORRELATION 


Figure 4. Comparison of computed and measured skin friction 
for an incompressible, axisymmetric , thin boundary 
layer in the absence of pressure gradient. 





4 . DISCUSSION 


The modifications described in Sections 2 and 3 have been made to 
the spinning-body version of the EDDYBL computer program. Con- 
sequently, the program now has provision for computing boundary- 
layer development segmented spinning bodies using the following 
turbulence models: 

1. Aguilar mixing-length model; 

2. Wllcox-Rubesin two-equation model; 

3. Wilcox-Rubesin RSE model. 

In the latter option either the original pressure-strain-correlation 
closure approximations or those devised by Launder, Reece and Rodi 
can be used. The program is now ready for general applications. 
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APPENDIX 

MODIFICATIONS TO INPUT AND OUTPUT 

In this Appendix we summarize all changes to program input and 
output which have been made in this research study. 

A.l INPUT MODIFICATIONS 

All modifications to input are in the input NAMELIST' s. The 
following input quantities have been added. 


NAMELIST SAFC0N 1 * 

VARIABLE DEFAULT 


NAME 

SYMBOL 

DEFINITION 

VALUE 

AHAT 

A 

a 

LRR Closure Coefficient 

0.76 

BHAT 

A 

6 

LRR Closure Coefficient 

0.11 

GHAT 

A 

Y 

LRR Closure Coefficient 

0.36 

RMIO 

‘Vi 

Radius of cylindrical body at which 
computation begins; to be used when 
C0NE <10~ 6 

1.00 

RSUBE 

R e 

Viscous Modification Coefficient 

1.00 

RSUBW 

R u 

Viscous Modification Coefficient 

2.67 

ZC1 

C 1 

LRR Closure Coefficient 

1.50 

ZC3 

°3 

Special Wilcox-Rubesin Closure 
Coefficient 

0.00 


^For clarity, the letter 0 is denoted by 0 in all F0RTRAN names 
while zero is denoted by 0. 
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NAMELIST NAM2 


VARIABLE 

NAME SYMBOL 


DEFAULT 

VALUE 

K£)DPRT 

Turbulent Prandtl number flag. 
K0DPRT = 1 gives Pr T * 0.89 
and K0DPRT * 2 gives 
Pr T - 0.95 - 0.*»5 (y/6) 2 

1 

KTC0D 

Transition width flag. 

KTC0D - 1 if TLNGTH input and 
KTC0D « 2 for TLNGTH computed 
from Equation (53). 

2 

SST s t 

Transition-point location.. 

l.E'8 

TLNGTH (s f /s t - 1) 

Transition width 

2. 


Using the default values gives the hybrid Wllcox-Rubesin and LRR mod 
el (WR#) described In Section 2. If the original WR model is 
desired, one uses the following values: 

a » 6 » 1/2, y = k/3 

R e - 1, R w = 2, C 1 = 9/2, C 3 = 5/2 

Note also that QBETA ($##) and SLAMDA (y) are no longer input but 
are computed according to formulae given in Section 2, viz, 

Equations (28) and (33). Also, note that in NAMELIST SAPC0N, the 
input variable NFLAG now has the following meaning: 


( ° 

, RSE Model 

NFLAG = J i 

, Two Equation Model 

I 2 

, Mixing Length Model 


Finally, note that for mixing-length runs (NFLAG =2), it is 
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necessary to impose a tighter convergence criterion than with 
advanced models. The following values are recommended: 

1 0.001 , NFLAG « 2 

0.010 , NFLAG = 0 or 1 

A. 2 OUTPUT MODIFICATIONS 

The only output modification is to the printed output which has 
been added for the mixing-length model. The printout at each 
station of skin friction, momentum thickness, etc., is as with the 
other two models with the following alternations: 

1. The quantities DSMX0, JPRINT and KPRINT 
are always zero; 

2. The quantity ES is the streamwise inter- 
mittancy, r(s). 

The first two sets of printed profiles are slightly different from 
those given for the advanced turbulence models. First the program 
prints n, y/<$, u/U e , y + , u + , (T/T e -1), T t /T t , (1 + e/v). Then, 
the program prints n» y/<5, V, w/U g , <-pu *v f >/p", c-pv'w'Vp, 
q/C_T U and 1/6 where A denotes mixing length and V is the 

P “ © 

2i 

dimensionless transformed vertical velocity. 
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